Critical conditions for a stable molecular structure 
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Here, we show how the molecular structure appears and becomes stable for supercritical physical 
conditions. In particular we consider, for the ammonia molecule in a gas, a simple model based 
on a standard non-linear double-well Schrodinger equation with a dissipative term and a term 
representing weak collisions. 

PACS numbers: 03.65.-w, 73.40.Gk 
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The existence of a well defined molecular structure for 
a symmetric molecule is an old and ongoing problem in 
chemistry ||^, Various attempts have been made to 
make quantum mechanics agree with apparently station- 
ary asymmetric states. It was clear from the beginning 
that the action of the environment would be the basic rea- 
son of this effect, since from a quantum mechanics point 
of view an isolated symmetric molecule has no structure. 
The point is to find the simplest model simulating the 
influence of the environment, and to exploit it. 

In 1927 Hund understood that in the case of cer- 
tain molecules, such as ammonia NH3, the basic model 
is a double well one, this simple model is able to give 
the splitting and the inversion frequency of an isolated 
molecule, but not the structure related to localization. 
From the experimental point of view, experimental data 
showing the decrease of the splitting, as a decrease 
of the inversion line (red shift), for increasing values of 
the gas pressure appeared in 1948. Theoretical mod- 
els for inter-molecular interactions were considered, and 
they show an interesting fit of the data at small pressure 

i, §■ 

Further researchers proposed some simple models in 
order to explain the localization phenomenon of a sym- 
metric molecule 0, [||. Later, some works ||^, [0 have 
exploited the strong effect of metastability, both static 
and dynamic, on red shift and localization. 

In 1995 it was shown |ll| that a non-linear 
Schrodinger equation, previously suggested by Pratt Q 
for the hydrogen case, is able to give spontaneous sym- 
metry breaking and the bifurcation of the ground state 
at the value 



1 



(1) 



of the dimensionless non-linearity parameter /i defined 
below (which is expected to be monotically dependent 
on the gas pressure). Later, it became clear that this 
non-linear model is able to give the red shift and local- 
ization too. In particular, in and Q it is suggested 
that localization appears for fi — 2/ic; while in p4[ it is 
suggested that the critical value (|^) found in |l^ coin- 
cides with the vanishing of the inversion line, and corre- 
sponds to the pressure of about 1,7 atmospheres at room 



temperature [ [l5| . 

In this letter we consider the non-linear double-well 
model as above [|l^ in order to obtain the stability of the 
localization for finite values of the parameters. Let us 
recall jl^ that actually we consider a substituted ammo- 
nia molecule NHDT, so that a localized state has one of 
the two possible chiral configurations. In the approxima- 
tion of a two- level system, an invariant quantity (energy) 
appears, and we observe periodic motions of two kinds: 
vibrational periodic motions around just one well (inside 
the grey regions of Fig. and ^) or beating periodic mo- 
tions between the two wells. In order to have a more 
physical model, we add a dissipative term, related to the 
photonic emission in the radio frequency range, giving the 
decreasing of the energy. We also take into account the 
effect of weak collisions. As a result we find that a chiral 
configuration is unstable if the non-linearity parameter 
is less than 3/ic- For larger values of this parameter, the 
chirality states become stable provided that the interval 
between two collisions is larger than the relaxation time, 
in order that when a collision occurs, the state be near 
an asymmetric stationary state. 

In order to better understand the situation, consider 
Figs. and ^ where the space of states is shown, as a 
sphere projected on a square by a Mercator map. The 
two coordinates z and 9 respectively represent the im- 
balance variable (which measures the localization) and 
the relative phase of the components of the state with 
respect to the two localized states. A collision actually 
changes 9 but not z. Thus, for fi greater than (see 
Fig. 1^) , the line defined by a value of z equal to the one of 
an asymmetric stationary state is fully contained in the 
vibrational region. In contrast, for /i < 3/ic (see Fig. ^) 
we have that a change of 9, due to a collision, could shift 
an asymmetrical stationary state into the beating region. 

The Hamiltonian for a single ammonia molecule takes 
the form Hq = —-^A + V where I^ is a double- well 
potential invariant under a coordinate refiection 

V = PVP, 

where P is the unitary symmetry operator representing 
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the inversion of the n-th coordinate 



(2^1, 



l), n>l. 



The tunneHng time through the inter-weU potential bar- 
rier is inversely proportional to the energy splitting AE 
between the odd- and even-parity eigenstates |— ) and 
1+) with energies A„ and A+, and the solution of the 
unperturbed equation 

zh^^m)) = Ho\m) 

shows a beating motion between the two wells with pe- 
riod 

The actual semi-classical parameter is the energy split- 
ting AE, and we choose the units such that h = 1 and 
AE < 1. 

Let us consider the interaction of a single molecule 
with the other molecules of the gas. In the mean field 
approximation, the new Hamiltonian takes the form 



H = Hn 



W 



where the term W is given by the polarization of the ex- 
ternal environment due to the presence of the single am- 
monia molecule itself. Because of the dissipative terms, 
the perturbation W could be written by means of a non- 
linear and non-Hermitian term 

W = jyg {el + ir^P) , = (VWIslV-W), e, < 0, 

where g{x) is a given odd function, PgP — —g, and where 
the parameters e and 77 respectively measure the strength 
of the dipole-dipole and the dissipation interactions. As 
in the case of complex Ginzburg-Landau equations fl^ 
and of Gross-Pitaevskii equations with weakly dissipative 
effects [|l^ , we expect to observe vortex solutions here. 

We underline that in the non-dissipative case, where 
77 = 0, then we have the conservation of the charge: 



m)\m)^\\mr = i 

and the conservation of the energy functional 



(2) 



£{ij)^\\Vil^f + {iP\V\^) + ^e{ 
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In the dissipative case, where t] < 0, since the term W 
is such that [TP, W] = 0, where T|V'(t)) = \^{t))* is the 
time-reversed operator, we still have the conservation of 
the charge (^) , and we expect that the energy functional 
decreases. The system then relaxes towards local minima 
of the functional. 

Here, we make use of the two-level model where we 
approximate the wave- function \tp{t)) by means of its pro- 
jection on the two lowest states. In particular, in such 



an approximation the total wavefunction of the system 
may be expanded as 



\m)^aR{t)\R)+ciL{tm, 
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where 



l^) = ^(l+> + |-», |i> = ^ (!+)-!-)) 

are the right and left hand-side states; they are such that 
P\R) — \L). The normalization condition (||) on the 
wavefunction implies that la^p + I^Lp = 1- 

By substituting \'^(t')) by (||) in the time-dependent 
Schrodinger equation 

it follows that the expansion coefficients a^j and have 
to satisfy to the following system of ordinary differential 
equations 

iixB. = ^a^. — loaL ^ evcaB.^ i"(]v^^L 
iiiL = ^a-L — ^o,R — evcaL — ir]iycaji (4) 

where 

c. = i(A_-A+), f] = i(A_ + A+) 

and where we set 

c^{R\g\R) = -{L\g\L). 

In order to study the beating motion, it is convenient 
to introduce the relative phase 

= arg(aij) - arg(aL), 

which is a torus variable, and the imbalance variable 

which takes values in the interval [—1,-1-1]. They have 
to satisfy to the system of ordinary differential equation 



z ~ ujZ{z,9) 

e = uje{z,e) 



(5) 



where 



Z{z, 9) = 2y/l - z2 [sin 6 ~ (^z cos 9] 



and 



Q{z,9) 



[cos 9 + (z sin 9] + 2fiz 



and 



/i = —c^e/uj, C — —c^rj/oj. 
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FIG. 1: In this figure we plot the parametric graph of the 
solution z{t), 6{t) of equation (5) for a given initial condition 
(point A) in the case = 1.5, C, = 0.2, point B denotes the 
position of the state &t t = lOr. Circle points denote the 
stationary solutions. In the non-dissipative case ((" = 0) 
we have only vibrational motions around one well inside the 
grey region, outside the grey region; we have periodic beating 
motions between the two wells. Here, z denotes the imbalance 
function taking values in the interval [—1, +1] and is a torus 
variable taking value in the torus (— tt, +7r]. 



Here, /i and C represent positive dimensionlcss parame- 
ters that measure the efFective non-linearity and the dis- 
sipation, respectively. 

We observe the existence of the critical value for 
the non-linearity parameter /i. We consider, at first, the 
weak non-linearity case such that ^ < Hc- In such a 
case, equation admits just two stationary solutions 
corresponding to the unperturbed even- and odd-parity 
eigenstates: zi = and 9i — 0, corresponding to the 
even-parity eigenstate, is a stable stationary solution, 
and 22 = and 62 — tt, corresponding to the odd-parity 
eigenstate, is an unstable stationary solution. 

Then, we consider the strong non-linearity case such 
that ^ > fic- We observe that when fj, takes the value fie 
then the stable stationary solution makes experience of a 
bifurcation phenomenon [ pl . More precisely, for fj. > fic 
we have 4 stationary solution; two of them still corre- 
spond to the unperturbed even- and odd-parity eigen- 
state, the other two correspond to asymmetrical states 
that, in the limit of large non-linearity, are fully local- 
ized on just one of the two wells. The stationary so- 
lution zi = and 9i = 0, corresponding to the even- 
parity eigenstate, is a saddle point for fj, > /ic; the sta- 
tionary solution Z2 = and 62 = t^, corresponding to 
the odd-parity eigenstate, is still an unstable solution; 
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FIG. 2: In this figure we plot the parametric graph of the 
solution for a given initial condition (cross point) in the case 
/i = 4.5, C, = 0.2. Diamond points denote the stationary 
solutions. 



2:3 — ^ (/i^ — 1)/ (/x^ -I- C^) and Bz — arctan(C23), and 
Z4 = — Z3 and 64 = —63 are stable asymmetrical station- 
ary solutions. Thus any state generically goes near to 
one of these two asymmetrical stationary states. 

Therefore, in the strong non-linearity case we have 
that any initial state, except the two even- and odd- 
parity unperturbed eigenstates, finally goes to one of 
the two asymmetric stationary eigenstates giving a chiral 
configuration for the ammonia molecule. 

Now, we show that this chiral configuration is sta- 
ble with respect to collisions when the pressure is large 
enough. We underline that in the ammonia case the 
thermal energy at room temperature is smaller than the 
distance between the doublet {X±} and the other energy 
levels, so that the validity of the two-level approximation 
holds, and it is much larger than the splitting energy, so 
that a collision could produce a strong variation of the 
energy Siip)- 

To this end we introduce a simplified model for molec- 
ular collision. When the single molecule undergoes a col- 
lision we add to the Hamiltonian H a perturbative term 
of the type f{x)v(t) where f{x) is a function with com- 
pact support and v{t) is a given time dependent function. 
For instance, let v{t) = XtiM [t) where % is the character- 
istic function on the interval [ti, t2\, and where we assume 
that the perturbation acts for a time much shorter than 
the beating period, that is 

^2 - ii < T, (6) 
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Equation then takes the form 



iciL 



since 



flafi — ujul + evcaj^ + irjvcaL + v(t)cRaii 
flaL — ujan — evcaL — irjvcaR + v(t)cLaL 
CR = {R\f\R), CL = {L\f\L) 



{R\f\L) = {L\f\R)^0 
Hence, the system (||) takes the form 
z = u!Z{z,6) 

9 = u;e{z,e) + {cR~CL)v{t) 

from which it foUows that \z{t2) ~ z{ti)\ ^ 1, since (|^), 
and where it is not possible obtain a similar bound for 
when cr^ CL- 

This fact, that is the relative phase is strongly mod- 
ified after a generic collision (such that cr ^ cl), does 
not actually destroy the chiral configuration of a local- 
ized ammonia molecule if the pressure is large enough, 
i.e. such that /i > S/ic- Indeed, in such a case we have 
that the stable solution (zs, ^3) (respectively {24, 64)) has 
a basin of attraction containing the strip z > z* (re- 
spectively z < —z*) if C is small enough and where 
z* = 2-\/2/3 < Z3. We explain this fact by means of 
a continuity argument in the limit of C = 0. Indeed, in 
such a limit we have the existence of two separatrix lines 
I p^ , [ p^ starting from the stationary solution (zi, 9i) and 
satisfying the equation 
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It is not hard to see (Fig. ||) that these paths are con- 
tained in the strip —z* < z < z* if /i > S/ic- As a result, 
it follows that a perturbation due to a collision acting in 
an interval of the order (^) shifts a state initially near to 
one of the asymmetric stationary stable eigenstates, to 
a state belonging to the basin of attraction of the sta- 
tionary eigenstate itself. In particular, the state will be 
always far enough from the unstable stationary state and 
thus, after a finite time depending on it returns near 
to the initial stationary state without visiting the other 
well, provided that in this period another collision does 
not occur. 

Finally, let us notice that the molecular structure is 
not completely stable for strong collisions; indeed we have 
the phenomenon of racemization which makes the statis- 
tical mean of the chirality vanish at large time, and that 
an external electromagnetic field in the radio frequency 
range could be able to destroy the molecular structure. 

In conclusion, in this letter we have shown that a dis- 
sipative non-linear model is able to explain the molecular 
structure of symmetric molecules. Such chiral configu- 
rations are stable for weak collisions provided that the 



non-linearity parameter is larger than the critical value 
3/j,c, and that the frequency collisions is small enough. 
The kind of trap we propose here is simple, but not triv- 
ial, and in any case is able to make stable the sponta- 
neous symmetry breaking given by the non-linearity. It 
could be relevant in the theory of dechoerence, related to 
the appearance of classical mechanics, and in the study 
of many irreversible precesses. As it clearly appears 
in this letter, we have only considered the spontaneous 
symmetry breaking effect. In fact, in the case of organic 
molecules, one enantiomer may be dominant. This ef- 
fect could be explained by means of a small initial enan- 
tiomeric excess, due to the parity violation for weak in- 
teractions , largely amplified during a very long time 
by the dynamics of the system j2^. 
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